DMRG study of the Bond Alternating S=l/2 Heisenberg ladder with 
Ferro-Antiferromagnetic couplings 



J. Almeida*, M.A. Martin-Delgado* and G. Sierra* 
* Departamento de Fisica Teorica I, Universidad Complutense. 28040 Madrid, Spain. 
'Institute/ de Fisica Teorica, C.S.I.C.- U.A.M., Madrid, Spain. 

We obtain the phase diagram in the parameter space (J' /J, 7) and an accurate estimate of the 
critical line separating the different phases. We show several measuments of the magnetization, 
dimerization, nearest neighbours correlation, and density of energy in the different zones of the 
phase diagram, as well as a measurement of the string order parameter proposed as the non vanishing 
phase order parameter characterizing Haldane phases. All these results will be compared in the limit 
J'/ J ^> 1 with the behaviour of the S = 1 Bond Alternated Heisenberg Chain (BAHC). The analysis 
of our data supports the existence of a dimer phase separated by a critical line from a Haldane one, 
which has exactly the same nature as the Haldane phase in the S = 1 BAHC. 
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I. INTRODUCTION 

Quantum systems when placed in low dimensional lat- 
tices typically exhibit strongly correlated effects driving 
them towards regimes with no classical analog. Many 
properties of these regimes or quantum phases pj de- 
pend in turn on the properties of their ground state and 
low lying energy excitations Q ■ 

A problem of particular interest in the field of strongly 
correlated systems is the emergence of critical phases in 
a system where the generic behaviour as coupling con- 
stants are varied is to be a gapped system, although those 
gapped phases may be of different nature. In this paper 
we address this problem by selecting a system of quantum 
spins that allows us to perform a detailed study of criti- 
cal and non-critical phases on equal footing, i.e., without 
any bias towards an a priori preferred phase. For reasons 
explained in Sect|IIJ the quantum spins are arranged in 
a 2-leg ladder lattice Q with anti-ferromagnetic Heisen- 
berg couplings along the legs while rung couplings are 
ferromagnetic. In addition, we also introduce an explicit 
dimerization coupling in the Hamiltonian along the leg 
directions, which can be varied from zero to strong val- 
ues. This coupling plays a major role in order to create 
the aforementioned critical phases out of a system with 
only gapped phases. 

This particular type of 2-leg ladder system has a num- 
ber of open problems such as the precise location of crit- 
ical phases in the phase diagram of the coupling con- 
stants, and the nature of the gapped phases it exhibits. 
Our study is complete enough so as to be able to solve 
for these problems in a very precise manner. 

The understanding of these purely quantum effects is 
usually a hard problem. Pcrturbative and variational 
methods in quasi-one dimensional systems like chains and 
ladders are not well suited to uncover the physics in the 
whole range of coupling constants involved in the descrip- 
tion of the interactions in the system. On the contrary 
the DMRG method Q, 0], @, 0, [| allows us to iden- 
tify the critical phases clearly and without any bias. This 



is so because the method is non-pcrturbative and allows 
a controllable management of errors. 



Our studies are also of interest since experiments on 
ladder materials have revealed a very complex behaviour, 
such as an interplay between a spin-gapped normal states 
and superconductivity 0]. Moreover, a new field of study 
for these complex effects has been opened by the simula- 
tion of strongly correlated systems in optical lattices [lj| , 
in particular quantum spin chains and ladders [ll[. 



This paper is organized as follows: in Sect|TT]we intro- 
duce the model Hamiltonian {1} describing a 2-leg ladder 
lattice of spins S = h with columnar bond-alternating 
antifcrromagnctic couplings in the horizontal direction 
and ferromagnetic couplings in the vertical direction, see 
Figjl] We can identify some particular behaviours in ap- 
propriate weak and strong coupling limits, but not for 
generic values of the couplings. In Sect Mil we point out 
the rich physical effects posed by open boundary condi- 
tions in these 2-lcg ladders with finite length, although 
it also implies an a priori analysis in order to find out 
which low-lying states contribute to the gap of the sys- 
tem in the thermodynamic limit. This we can be done 
with the DMRG method by targeting several states and 
measuring their magnetization properties in the bulk and 
at the ends. Then, we compute numerically the gap and 
we establish the existence of a critical line in the quan- 
tum phase diagram of the model. A numerical fit of this 
critical curve is also given. In Scct llVI we determine the 
structure of the phase diagram by identifying the type of 
gapped phases occurring at each side of the critical line 
found in the previous section. They correspond to Hal- 
dane and dimer phases. They arc identified by measuring 
the string order parameter and the dimerization param- 
eter with the DMRG method. We complete our study of 
these phases measuring different obscrvables. Sect|V] is 
devoted to conclusions. 
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FIG. 1: Pictorial representation of the quantum Hamiltonian 
([T]). The geometry of the lattice is a 2-leg ladder. Each solid 
dot is a spin S — 5. In the horizontal direction (legs), we 
picture the bond alternation with strong links J(l +7) and 
weak links J(l — 7). In the vertical direction (rungs), the 
system is arranged in the form of a columnar dimerization: 
strong links are parallel to one another, and similarly for weak 
links in the lattice. 



II. THE MODEL 

Competing ferromagnetic versus antifcrromagnctic 
spin interactions may give rise to critical phases if they 
are appropriately arranged in certain quasi-one dimen- 
sional lattices. One emblematic example of this phe- 
nomenon is a lattice of quantum spins with the shape of 
a 2-leg ladder such that there are antiferromagnetic cou- 
plings along the legs and ferromagnetic interactions along 
the rungs connecting both legs. In addition, the antifer- 
romagnetic couplings are bond-alternating in a columnar 
fashion. Dimerization interactions in the Hamiltonian 
are also known as staggered interactions. This configu- 
ration is shown in Fig[TJ More precisely, this configura- 
tion of Heisenberg-like interactions is associated with the 
following quantum Hamiltonian 



L-l 



£=1,2 i=l 
L 



(1) 



J'5^S<(1)-S 1 (2), 



i=l 



where Si(£) are quantum spin 5 = | operators located 
at site i of the leg £, and J > 0, J' < 0, 7 E [-1, 1] 
are the antifcrromagnctic, ferromagnetic and staggering 
couplings, respectively, as mentioned above. 

Notice that several known regimes can be reached by 
tuning the coupling constants towards particular val- 
ues. In the weak coupling limit, making \J'/J\ -C 1 
we end up with a system consisting on two effectively 
decoupled 5 = 1/2 Hcisenbcrg chains with bond alterna- 
tion (BAHC), which are known to be gapped for every 
value of the dimerization parameter 7 [12j, except for 
the point 7 = 0. In the strong coupling limit, making 
\J'/J\ 1, J' < the system can be effectively de- 
scribed by a 5 = 1 spin chain with bond alternation, 
which is predicted to be gapped for all values of 7 except 
for a critical point at a non-zero value 7 C fl3j | . These pre- 
dictions are based on an approximate mapping onto the 
0(3) a model Q at topological angle 9 = 2itS(l — 7). 



similarly another symmetric critical value at j c = — 5. 
Thus, we shall always concentrate in the region 7 > 0, 
due to the symmetry 7 <-> —7 in the Hamiltonian (JTJ) . 
This non-linear sigma model (NLctM) prediction misses 
the correct location of the critical point due to the ap- 
proximations involved in that mapping. The exact loca- 
tion of this point has been widely studied [l6| and results 
slightly varied depending on the approach, however mod- 
ern studies place it at j c = 0.259 [la] ' a l so compatible 
with Figg]>i/;er) which gives j c = 0.2590 ± 0.0001 for a 
chain of 500 sites. These studies also conclude that the 
region | — y | < 7 C corresponds to a Haldane phase while for 
7 1 > 7 C we move to a dimer phase. The emergence of a 
dimerized 5 = 1 spin chain in the strong coupling limit 
can be explained by noting that as the rung coupling is 
ferromagnetic and strong J' < 0,\J'\ J, the two spins 
5 = I in each rung find energetically favorable to form 
a spin triplet. 




8.0 



perpendicular Hm 

odd x ■** x 
cumulative sum a x x * *. x x ^ 



-0.4 




-0.4 



position in the chain 



This yields a critical value of 7 C 



when 9 — it. and 



FIG. 2: Computations on a L = 2 x 96 ladder at the point 
( J 1 1 J = —2.5, 7 = 0). Up: Mean magnetization (Sf) in the 
states with St ot = 1 and Sf at — 2, as explained in the text. 
Inset: Cumulative sum of the magnetization over the whole 
extent of the ladder and the states with Sf ot = 1 and Stot — 2. 
The order of the sites in the a;-axis corresponds in this case to 
the path used to traverse the ladder in a DMRG sweep. Down: 
Difference of energy density (DOE) of the excited states with 
5 t z ot = 1 and S t z ot = 2: (S/,iS// ii /)5*=a - (Sf, j S^ ji /)s^=i- -The 
scale on the right axis corresponds to the cumulative sum. See 
text for more explanations. 

For generic values of the coupling constants in the 
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FIG. 3: Gap A21 computed on a L = 2 x 140 ladder for 
different values of the parameter J' / J. Each minimum in the 
gap belongs to the critical line. 



and only in the thermodynamic limit L — > oo it may 
vanish for certain values of J' / J and 7 which define the 
critical line we are searching for. Thus, the signature 
of a gap in A(J'/J, 7) for fixed J'/ J and varying 7 will 
show up as a minimum in the dimerization parameter. 
Upon increasing the value of L, we shall obtain more 
robust estimations of the critical value 7 C (J'/J) from the 
minima 7 m i n (X). This is a finite-size scaling analysis of 
the DMRG numerical data. 







Hamiltonian (JTJ) , this model has been the subject of a se- 
ries of conjectures based on exact diagonalization numer- 
ical studies [H| in the absence of dimerization 7 = and 
analytical studies using bosonization and NLctM map- 
ping [il^ ] in the presence of dimerization 7^0. Those 
numerical methods only allowed to reach ladder lengths 
typically of L = 6 or so, which prevents from reaching 
any definitive conclusion on the bulk properties of the 
system in the thermodynamic limit. As for the analyti- 
cal studies, they conjectured the existence of a possible 
critical region, but due to the nature of the methods it is 
not possible to give its location in terms of the original 
coupling constants in the model Hamiltonian (fTj). 



III. CRITICAL REGION 

One of the main issues in this model ((T|) is whether it 
exhibits a critical line in the quantum phase diagram of 
J' / J vs. 7. We solve this open problem in the positive 
by using the DMRG method in finite version algorithm 
which provides us with better accuracy values than the 
infinite method version, although at the expense of more 
demanding time computing resources. The performance 
of the finite DMRG algorithm is characterized by the 
following parameters: the number of states m retained 
in the truncation process of the RG method, the weight 
of the discarded states w m which is a measure of the 
DMRG error, the number of sweeps n s or iterations of 
the method after the initial warm-up process and the 
tolerance e of the target state energy which controls the 
average number of iterations that will need the diagonal- 
ization algorithm (Lanczos in our case) to compute the 
target state. We shall provide values of these parameters 
in our numerical computations below. 

Before applying the finite-size DMRG method, two im- 
portant remarks are in order: 

i/ As we shall always work with a fixed value of L the 
length of the lattice, the gap A(J'/J, 7) is always finite 
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FIG. 4: Tbp:Critical region of a ladder of size L = 2x 140. 
The points of the critical line correspond to the coordinates 
that minimize A2i(7, J' /J). The solid line is only a guide for 
the eye. Bottom: The value of the minimum gap for the ladder 
with L = 2 x 140 sites is very similar to the corresponding 
S — 1 BAHC with L = 140 but the value of y c that minimizes 
this gap is still a bit shifted, which constitutes a signal that 
J' /J = — 20 is still a low value to accurately mimic the limit 
BAHC behaviour. The computations for the L = 500 BAHC 
were performed storing m — 450 eigenvectors of the density 
matrix. 

ii/ The physics of this 2-leg ladder (p} is richer when the 
lattice has open boundary conditions. Moreover, the nu- 
merical performance of the DMRG method is also better 
under these conditions. However, open boundary condi- 
tions must be handled with care in order to identify the 
gap A(J'/J, 7) we are after. We shall provide ways to 
do this identification by targeting appropriate low-lying 
states and measuring convenient observables with them. 

In particular, using open boundary conditions we have 
found that the first excited state lies within the sector 



4 



with total z-spin angular momentum S* ot = 1 , but in the 
Haldane phase it converges to the ground state that has 
5f ot = as we take larger sizes of the system. We have 
then to consider the next excited state in the sector with 



2 to compute the gap of the spectrum as 
A 21 ~E (S^ t = 2)-E (S^ t = l), 



(2) 
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FIG. 5: The region of the critical line in the limit \J'/J\ 3S> 
1 fits very well to a potential function of the form J'/J = 
(7(0.259 - 7 c ) fe , with C = -0.16 ± 0.01 and k = -1.25 ± 0.01. 

The reason for considering A21 instead of A10 : — 
E (St ot = 1) — Ea(S* ot = 0) as the gap of the system, can 
be justified as follows: in the complete dimerized limit 
7 = 1 , it is clear that the difference in energy between two 
arbitrary consecutive levels is the same, and corresponds 
exactly to the energy needed to promote one singlet bond 
to a triplet. The argument for the limit 7 = makes use 
of the properties of Haldane phases, where it is known 
to appear a non bulk excitation due to the existence of 
virtual spins at the end of the chain. Our conclusion 
is that the lowest lying state with S"* * = 2 consists on 
a superposition of two kind of excitations, namely the 
Haldane non-bulk triplet mentioned before and the bulk 
itself, also giving a triplet. Considering this scheme, in 
order to obtain the gap related to the bulk excitations 
we have to substract the non bulk excitations present in 
the lowest lying states of sectors Sl ot = 1 and 5' ot = 2. 

In Fig|2] we show rigourous comparations of these 
two states in the Haldane limit 7 = 0. Computations 
have been done on a L = 2 x 96 ladder at the point 
(J'/J = —2.5,7 = 0). On the up part of the figure, 
we plot the mean magnetization (S?) in the states with 
Sf ot = 1 and Sf ot = 2, computed in one leg of the 
ladder, since due to the symmetry of the Hamiltonian, 
the magnetization is the same in both legs. As a check 
of the accuracy of our computations we observed that 
the results in both legs are the same up to the fifth or 
sixth decimal digit. In the Inset of that figure, we show 
the cumulative sum of the magnetization over the whole 
length of the ladder and the states with S* ot = 1 and 
S* ot = 2. The order of the sites in the x-axis corresponds 
in this case to the path used to traverse the ladder in a 
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FIG. 6: String order parameter computed for a ladder con- 
sisting on L = 2 x 96 sites. The SOP has been computed 
forming the triplets with adjacent S — 1/2 spins located in 
different legs. Clearly the value of this parameter is non- 
vanishing in the low region of 7, where the system is in the 
Haldane phase. The dimer phase is nonetheless character- 
ized by a vanishing SOP Jnset: SOP computed for a S = 1 
BAHC. The resemblance between both systems is evident in 
the region with \J'/J\ > 1. 



DMRG sweep. In the down part of this figure, we plot 
the difference of energy density of the excited states with 

^tot = 1 an d S'tot = 2: (Se,iSe',i')s*=2 — (S£,iSi',i')s*=i 
The difference has been divided into three contributions: 
the contribution labelled with even stands for links in- 
volving sites in the same leg and the even sublattice 
(£' = l. i = 2k, i' = 2k + 1), odd involves links joining 
sites in the same leg and the odd sublattice (£' = £,i = 
2k — l,i' = 2k), and perpendicular denotes links among 
legs (£ = 1,1' = 2, i = k,i' = k). The cumulative sum 
of the difference of the various contribution, measured in 
the right axis scale, is also shown. Interestingly enough, 
we can observe the magnetization pattern at the ends 
being almost identical in the states with S% ot = 1 and 



2. The contribution to the z-axis projection of 



the spin coming from the ends is equal to 1 in both cases. 
Notice also that the difference of the density of energy be- 
tween these states is close to zero at the ends, while it 
becomes clearly apreciable in the bulk. All these facts 
strongly support the picture of a non-bulk triplet exci- 
tation with the same nature in both states, that leaves 
the bulk of the chain with a neat value of the projec- 
tion equal to S^ ulk = and S^ ulk — 1 and gives a strong 
hint on the equivalence of A 10 with periodic boundary 
conditions and A21 for open systems. 

After this previous analysis to identify the states 
needed to target the gap of the system, we present in 
FigE] some values of the gap A21 for a ladder consisting 
on L = 2 x 140 sites, at different regions of the parameter 
space. Computations have been performed retaining 
m = 300 states of the density matrix and the grid used 
to explore the phase diagram is 7 G [0, 0.4], and — J'/J = 
{0.00, 0.20, 0.40, 0.60, 1.25, 2.50, 3.00, 3.50, 4.00, 4.50, 5.00, 
6.00, 7.00, 10.00, 15.00, 20.00}. The existence of a set 
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FIG. 7: A picture of the path used to compute the string order 
parameter in the 2-leg ladder with columnar bond alternation. 
The ellipses mean that the sites within them are forming a 
triplet. 
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FIG. 8: Different values of the dimerization parameter 
Di X^=i,2( s ^-i s <^) ~~ (Sf,iSf, i+ i) computed in the mid- 
dle of a ladder of L = 2 x 96 sites. Since we are explic- 
itly introducing some staggering in the Hamiltonian ([TJ), the 
dimerization parameter is non- vanishing even in the Haldane 
phase. However, the shape of the graphs seem to have an in- 
flexion at the critical point. The graph corresponding to the 
5 = 1 BAHC, has been scaled down by a factor 1/2 due to 
the effective coupling constant of the ladder, which is known 
to be half the constant corresponding to the BAHC. 



of minima in the function A2i(7, J'/J) is clear in this 
graph, although they shall become more distinguishable 
as we move to higher values of \J'/J\. 

As an instance of the accuracy of our results, we point 
out that a systematic examination of the error in each 
of the truncations of our DMRG computations reveals 
that the highest values in the whole process are of the 
order of w m ~ 10 -8 , and mostly they are of order 
w m ~ 1CP 10 . To obtain a suitable acuracy in the re- 
sults we have set the number of sweeps n s = 2 and the 
tolerance to le — 9. To compute with enough precision 
the critical value 7 C (J'/J) that minimizes the gap it be- 
comes necessary to use large amount of data. On this 
regard, we have used interpolated values resulting from 
the DMRG computations. 

Now, we can detect the presence of a critical line in 
the quantum phase diagram separating gapped phases. 
In FigfJ] we plot the critical region consisting of the co- 
ordinates for each minimum in Fig[3J In earlier studies 
we placed the critical point of the 5 = 1 bond- 
alternating Hciscnberg chain (BAHC) at j c = 0.259. The 
curve shown in Figd] shows a vertical asymptota that is 
still a bit off from this limiting value corresponding to 



the region \J'/J\ » 1, but this is simply because we 
have chosen a value of J'/J = —20 which is still not big 
enough and also due to finite-size effects on the 2-leg lad- 
der. In the lower plots of Figf5]we address these possibil- 
ities by comparing our ladder in the strong ferromagnetic 
limit with a pure 5 = 1 BAHCs with different sizes. Two 
parameters are important in this discussion, namely, the 
value 7c(J'/J) that minimizes the gap, and the value of 
the gap itself at this point A2i(7 c ). As we can observe in 
FigEl the value of A 2 i(7 c ) does not strongly depend on 
the particular choice of the coupling constant ratio J'/J, 
while it is definitely influenced by the size of the system. 
In FigUlower) it is shown that the shift of j c computed 
for two 5=1 BAHCs with different sizes, but still large 
enough both, is less noticeable than the difference in their 
value of A2i(7 c ). It is clear the similarity of this mag- 
nitude in the case of the ladder and the corresponding 
BAHC, as well as the shift in the value of j c . All this 
make us conclude that in order to attain a better con- 
vergence with the 5=1 BAHC and a better estimate of 
the critical asymptota j c = 0.259, we shall increase the 
strength of the ferromagnetic coupling rather than the 
size of the system. 

As we have a set of numerical data from the finite- 
size analysis of the critical line, we can also make a nu- 
merical estimation of the criticality curve. In fig. [5] we 
present a fit of the critical curve in the region close to 
7c ~ 0.259. We choose as trial function for this fitting an 
inverse power law with some coefficients and exponents 
that are fixed by our numerics, namely, 

The fitting yields the following estimations for the values 
of the parameters C and k that best fit the data: C = 
-0.16 ± 0.01 and k = 1.25 ± 0.01, and for simplicity the 
value of 7 C = 0.259 is taken for fixed. 



IV. HALDANE AND DIMER PHASES 

Once we have established the existence of a critical 
line in the quantum phase diagram of the model (fTJ), it is 
natural to wonder about the two gapped phases that this 
line separates. More specifically, whether they are differ- 
ent or not and their identification as quantum phases in 
the framework of strongly correlated systems. 

The possible nature of those phases can be guessed 
from the strong ferromagnetic limit 1 <C | J'/J\ of the lad- 
der, effectively leading to the 5 = 1 BAHC. The phases 
of this chain are known to be the massive Haldane phase, 
separated by a critical point from the also massive dimer 
phase. To test the nature of each phase, we will resort 
to two different order parameters. The Haldane phase 
is known to exhibit a particular hidden order that can 
be measured by the string order parameter (SOP) [2fl |. 
[2l| . The definition of this operator for a spin-1 chain is 
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FIG. 9: Correlation (S- (1)S?(2)} in the ground state of a 
L = 2 x 96 sites ladder and J' /J = -20. The plot on the 
left shows the correlation value for some arbitrary values of 
7. The plot on the right shows the value of this magnitude in 
the whole region of the parameter 7. 



as follows: 



o(i) = {siYl e *" si si 

k=2 



(4) 



This operator acting on our ground state measures how 
far it is from a spin liquid Neel state consisting on a 
sequence of S = 1 spins such that every spin with pro- 
jection Sf = ±1 is followed by Sf +k = =Fl and Sf +k , = 
for every < k' < k. 

When we deal with S — 1/2 particles, to compute 
the SOP we have to define the pairs of particles which 
are most likely to couple to give a triplet and compute 
the SOP along the path connecting them. In our case, 
the existence of a ferromagnetic coupling clearly suggests 
that the triplets will result via this coupling. It is also 
worth recalling that the SOP is a parameter suited to 
work with translational invariant systems. In order to 
correctly estimate the SOP in open systems, we must re- 
strict the computation to a region shorter than the whole 
length of the chain where end-effects are negligible and 
only bulk physics is relevant. In FigJSjwc show the SOP 
computed traversing the path shown in Fig(7j We can ob- 
serve a non-vanishing SOP in the Haldanc region, while 
it rapidly decays to zero in the dimer phase. The inset 
shows the SOP computed for a S = 1 BAHC and the 
resemblance between both systems is apparent. 

Therefore, the phase below the critical line in the 
numerical phase diagram of Fig[5] is a gapped Haldane 
phase. 

As for the region above the critical line in FigHJ we 
have guessed from the strong coupling limit that it may 
be a dimer phase. The structure of a dimer phase is such 
that full translational invariancc symmetry of the system 
is broken by one unit cell of the lattice. This situation 
can be detected by means of the dimcrization parameter, 
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FIG. 10: Correlation (5? (1)5? (2) > in the first lying excited 
state in the sector with 5 t z ot = 1 of a I = 2 x 96 ladder and 
J'/J = -20. 
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FIG. 11: Correlation (Sf(l)Sf (2)) in the first lying excited 
state in the sector with Stot = 2 of a ladder consisting on 
L = 2 x 96 sites and J'/J = -20. 



which can be defined for our particular 2-lcg ladder as 
A := < S >-iW ' S *W> - < S *W • S .+iW)- ( 5 ) 

£=1,2 

The subindex i is necessary since open systems are intrin- 
sically not translationally invariant. In FigJS]it is plotted 
the dimerization parameter of the ladder measured in 
the middle of the chain. Since staggering is explicitly in- 
troduced into the Hamiltonian the order parameter 
vanishes only at 7 — 0, but is finite even in the Haldanc 
phase. Nevertheless, our plots clearly exhibit different 
behaviours related to the convexity of the parameter at 
each phase. This observation indicates that an accurate 
estimation of the point of inflexion in the dimerization pa- 
rameter could be used as a measure of the critical point 
separating both phases. 

We have also performed some measurements in the 
ladder to give more hints to understand the nature of 
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FIG. 12: Magnetization curves in a L = 2 x 96 ladder for the 
lowest lying states in the sectors 5f ot = 1 (top), and t3f ot — 2 
(bottom). The curves are separated in the different sublat- 
tices consisting on the sites ocuppying odd or even positions. 
Notice that both states present a peaked magnetization at the 
ends for 7 = 0.0 well into the Haldane phase, while it vanishes 
in the dimer phase with 7 = 0.4. 



both phases. Figures l9l ITOl and [TT] show the correlation 
(Sj(l)Sj(2)) between sites in the perpendicular rungs. 
The pattern of the correlation can be understood by 
noticing that the correlation between two isolated S = 
1/2 spins coupled to give a singlet is (SiS 2 }/3 = —1/4 
while it equals 1/12 if the spins form a triplet. From 
these values we observe that the perpendicular rungs in 
the ground state are forming triplets and the distribu- 
tion is uniform all along the ladder. In the excited states 
however, the triplet strength of some rungs is weakened, 
signaling the presence of magnon-like excitations, also 



apparent in Fig[12] The nature of the non-bulk exci- 
tation present in the Haldane phase is not magnon-like 
and that explains the different number of kinks in the 
Haldane and dimer phase in FigrTHl and FigfTTI 



V. CONCLUSIONS 

We have determined the existence of a critical line 
in the quantum phase diagram of a 2-leg ladder with 
columnar dimerization and ferromagnetic vs. antiferro- 
magnetic couplings. In this study, we use the finite-size 
system DMRG method which allows us to give the lo- 
cation of that critical curve. Moreover, we have clearly 
identified the two phases separated by the critical line 
to be a Haldane phase and a dimer phase. This iden- 
tification is carried out by measuring the string order 
parameter and the dimerization order parameter in the 
whole range of values of the coupling constant ratio J'/ J 
and dimerization parameter 7. 

As a byproduct, we have introduced a systematic anal- 
ysis of the spins at the borders of the open 2-leg lad- 
der lattice. Our model is based on S = | spins, then 
these end-chain spins exhibit physical effects of their own. 
They are real spins unlike the virtual spins appearing 
in integer spin chains or ladders. Their physics is spe- 
cially interesting when the system size is finite, and even 
during the process of reaching the thermodynamic limit 
they produce non-trivial finite-size effects along the way. 
These facts difficult the technical analysis of the open- 
ing or closing of a gap in the low-lying spectrum of a 
2-leg ladder with open boundary conditions. We have 
solved these difficulties by analyzing the ground state 
and low-lying energy excitations with respect to their 
bulk and boundary properties such as local magnetiza- 
tion and the like. With this information, it is possible 
to identify which states contribute to the gap in ther- 
modynamic limit. These low- lying states have a definite 
total spin S z and they can be targeted with the DMRG 
method. In this fashion, we have been able to identify 
the gapped or gapless behaviour of the model within the 
framework of the DMRG with open boundary conditions. 
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